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Abstract 

The complete one-loop self energies (real and imaginary parts) for photons, 
gluons, electrons and quarks at finite temperature are calculated numerically 
and compared to the results of the hard thermal loop (HTL) approximation 
used for the resummation technique of Braaten and Pisarski. In this way some 
light is shed on the validity of the weak coupling limit assumption {g 1) 
or equivalently the high temperature assumption, on which the HTL approx- 
imation is based. Furthermore, the gauge dependence of the fermion self 
energy beyond the HTL approximation is considered. Finally the dispersion 
relations following from the real part of the self energies are compared to the 
HTL results. 
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I. INTRODUCTION 



In relativistic plasmas, such as a quark-gluon plasma (QGP) possibly formed in rela- 
tivistic heavy ion collisions or a QED plasma in a supernova explosion, and in relativistic, 
degenerate matter, as e.g. strange quark stars, medium effects are of great importance. For 
example, the originally massless or light particles of the system acquire an effective mass 
caused by the interactions with the other particles of the system. Also screening effects due 
to the presence of (color) charges in the system (Debye screening) prevent or at least reduce 
infrared singularities coming from the long range interactions in QED or QCD. 

These effects can be calculated from the self energy of the particles. Evaluated perturba- 
tively to lowest order at high temperature T or density (corresponding to a large chemical 
potential fi) the effective mass and the Debye mass are of the order gT or gfi, where g is the 
coupling constant. These masses can be large compared to the bare masses, which are neg- 
ligible at high temperature. The dispersion relations of the quasiparticles thus deviate from 
the ones of bare particles leading to important modifications of the equation of state . 
Therefore medium effects are an essential feature of relativistic systems at finite temperature 
or density. 

So far these medium effects have mainly been studied by considering the self energies in 
the high temperature approximation or, equivalently, in the hard thermal loop (HTL) 
limit 1^. In this approximation the momenta of the internal particles of the self energy 
are assumed to be hard, i.e. of order of the temperature or chemical potential, whereas 
the external momenta are soft, i.e. of the order gT (gfi). These momentum scales can be 
distinguished only in the weak coupling limit (7 ^ 1, on which the HTL approximation 
is based. The reason for these restrictions is twofold. First this approximation allows for 
analytic expressions of the self energies. Secondly, which is more important, the self energies 
calculated in this way are gauge independent guaranteeing gauge independent results for 
the effective masses, the dispersion relations, and the Debye screening. As a matter of fact, 
the HTL approximation builds the basis for a consistent extension of the naive perturbation 
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theory, which suffers from infrared singularities and gauge dependent results for observables. 
Resumming the HTL self energies (and vertices) leads to effective Green's functions, which 
have to be used in the case of soft external momenta (Braaten-Pisarski method 0). In this 
way consistent results for physical quantities can be found, which are gauge independent 
and complete to leading order in the coupling constant. At the same time the infrared 
behavior is greatly improved due to the use of resummed propagators, although there might 
be still infrared divergences owing to the absence of a static magnetic screening in the HTL 
self energies of the gauge bosons. For applications of this powerful method to interesting 
problems of the QGP we refer to Ref. 0. 

Here we want to compute the self energies and dispersion relations of photons, gluons 
and fermions at finite temperature but zero chemical potential within the one-loop approx- 
imation abolishing the HTL assumption and compare the results to the ones of the HTL 
approximation. This will give us a hint for the validity of the weak coupling assumption, on 
which the HTL approximation and thus the Braaten-Pisarski method relies. Up to now there 
are only explicit expressions for the gluon self energy beyond the HTL limit, where the next 
term of the high temperature expansion has been evaluated For a complete one- loop 
calculation, which can be done only numerically, there exist only integral expressions 0, 
results for certain momentum limits [|lO|, or analytic expressions for zero temperature but 



non-zero chemical potential |jTT|. Furthermore, it is instructive to have a look at quantities 



which vanish at leading order of the one-loop high-T expansion. As examples we consider 
here the imaginary parts of the self energies above the light cone and the gauge dependence 
of the fermion self energy. 

II. PHOTON SELF ENERGY 

First we consider the photon self energy or polarization tensor of QED within the one- 
loop approximation. In the vacuum it reads {K = (fco,k), k = |k|) 

^?riP) = -^e' J ^tr[7,5(-Q)7.5(ir)] (1) 
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where S denotes the electron propagator and P = K+Q. In order to compute H^^^ at finite 
temperature we apply the imaginary time formalism and, for temperatures T ^ rrie, neglect 
the electron mass, i.e. S{K) = I^Ap{K) where Af{K) = l/K"^ with ko = i{2n + l)7rT for 
integers n. To perform the sum over the Matsubara frequencies /cq it is convenient to use 
the Saclay representation of the propagators |T^: 



dre^''^ Af{t, k) 



(2) 



with 



A^(r, k) = ^[[l- riFm^'^^ - nF{k)e^^} 



(3) 



where nF{k) = l/[exp(A;/T) + 1] denotes the Fermi distribution. 

Defining n^^^ = U^ff^ and U^^^ = \{Sij — piPj)Ilfj^^ , [p = p/p), the longitudinal 
and transverse projections, respectively, of the self energy read 



n?^^(P) = -46^ ^T|:(Mo + k-q)Ap(g)A^(ir) 

n«^^(P) = -Ae" I -0-^Tj2ikoqo - (p ■ k)(p ■ q)) A^(g)Ap(ir) . 



(4) 



In the HTL approximation, i.e. assuming pQ,p -C T, analytic expressions can be obtained 



nr^(P) 



1 _ In P^^P 
2p po-p 



2 ' p'^ 



1- 1-^U^^ln^ 
V Pt 2p Po-p 



(5) 



where = eT/3 denotes the effective, thermal photon mass. The derivation of (§) can be 
found, e.g., in the appendix of ref. 0. 

Now we want to evaluate the complete one- loop self energies starting from (1), which we 
will write as 



d^k 

d^k 

(2^ 



Ff,f^ ~I~ k ■ qP^^^T^ 

p(;].^)-(p-k)(p-q)p(°].; 
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(6) 



where 

ko 

= \\[l- npiq) - npik)] ( \ \ ) 

+ [nM-nFik)](^ \—-]] (7) 

\Po + k-q Po-k + qjj 

and 

ko 

= zM [1 - - ^^(^)] f \ ) 

4:kq [ \po + k + q Po - k - q J 

-[nM-nF{k)](^ T— II- (8) 

\Po + k-q Po-k + qJ } 

Next we neglect the vacuum parts, i.e. the terms which do not contain any distribution 
functions in (|^ and as we are interested only in the temperature dependent medium 
effects. The angle integration over = p ■ k can be performed analytically. (For this purpose 
we remove the ?7-dependence in npiq) by shifting k ^ p — k in the corresponding terms.) 
The remaining integral over k has to be done numerically. The final results for the real 
part of the self energy are shown in Fig. 1 and Fig. 2, where they are also compared to the 
HTL results (^). For the longitudinal component, we observe a surprisingly good agreement 
between the two approaches even when the HTL assumption pQ,p <^T is not fulfilled. The 
transverse component, however, is more sensitive to the violation of the HTL assumption. 
For small p and po the one-loop result approaches the HTL limit as expected. Obviously, 
near the light cone the results of both approaches coincide even outside the range of validity 
of the high-T assumption. Indeed, the behavior of, e.g., the longitudinal self energy, for 
Po ^ P and all values of p is readily found to be H^'^'^ ~ 3m^po/(2p) ln[(po + p)/{Po — p)], 
which is just the HTL result (P) extrapolated to arbitrary, nearly light-like momenta. 

Next we consider the dispersion relations, which are given by the poles of the effective 
photon propagator following from resumming the self energy. For the longitudinal mode 
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(plasmon), which is absent in the vacuum, and the transverse mode, respectively, they are 
determined by 

l)*^QED^pyl ^p2_ nQ^^(p) = . (9) 

As in the vicinity of the hght cone the deviations of both approaches are small, the same is 
expected for the dispersion relations for not too large couplings. Indeed, even in the case 
of extrapolating the coupling constant to e = 3, the difference is less than 5% (Fig. 3). We 
note that in a good approximation the one-loop dispersion relations follow /^from the HTL 
approximations by a suitable rescaling of m^. As expected, both branches meet for p = 
in either approximation. 

III. GLUON SELF ENERGY 

Here we compute the complete one-loop gluon self energy which, in contrast to the abelian 
boson polarization, is gauge dependent. In Feynman gauge, the pure gauge contribution of 
the gluon loop, the gluon tadpole and the ghost loop diagram is 

— — [pV - + 2{KQg,, - K,Q, - K^Q,)] Ab{K)Ab{Q) 

(10) 

where Ab{K) = l/K"^ with = i2mTT. The Saclay representation of the bosonic propaga- 
tor reads 

r/3 



ABiK) = - f dTe'°-AB{T,k), (11) 

JO 



with 



AB(r, A;) = ^ {[1 + usme''^ + ^B(fc)e'"} (12) 

and the Bose distribution nsik) = l/[exp(A;/T) — 1]. As for the photon, we consider the 
longitudinal and transverse components of U'^^'^. Expressed by the functions Fj^'^\ 
defined in analogy to the fermionic functions in (|^,§) with —np replaced by n^, they read 



nf" = ^9'l^s [<b' + (2k ■ q + 



(2vr) 

n?"^^ = / ['^S^' + - - 2(P • k)(p . q)) . (13) 

Under the HTL assumption, terms oc Po'P^ drop out and the integrals over the bosonic 
functions yield up to a factor of two the same expressions as for the fermionic case, so the 
self energies obviously approximate to just the abelian expressions @ with now replaced 
by the thermal gluon mass nig = gT / -\/3- 

In comparison to the abelian case for the same values of pcP, the difference to the HTL 
result is evidently larger (Fig. 4 and Fig. 5). Again we find that for soft momenta as well 
as for po ^ p the HTL limit is approached and that the transverse part is more sensitive to 
the overstraining of the HTL assumption. 

In general, beyond the leading contribution of the high-T expansion, the (one-loop) self 
energies have a non- vanishing imaginary part above the light cone, for po > p. For the 
gluon self energies it is shown in Figs. 6 and 7. While the general requirement Q'(H) < 0, 
corresponding to a positive damping rate, is fulfilled for the one-loop photon self energy, it is 
violated by, e.g., the transverse gluon self energy (in Feynman gauge). Q This inconsistency 
indicates that extending the high-T expansion beyond leading order requires higher-loop 
contributions to obtain a consistent propagator. The positive imaginary parts above the 
light cone thus provide a measure of the magnitude of these sub-leading terms which may 
be of a similar order as the imaginary HTL contributions for po < p. 

Although self energies containing gluon propagators depend on the choice of the gauge, 
the dispersion relations are, in principle, observable and thus gauge independent. While 
the gauge dependent terms of the gluon self energy are (somewhat surprisingly) found to 
be of sub-leading order in the high-T expansion and the HTL-dispersion relations are thus 
gauge invariant, this is not the case in the one-loop approximation. Postponing the question 



^ Historically, this was one of the main triggers for the development of the Braaten-Pisarski 
method. 
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of the gauge effects to the next section, here in Feynman gauge the influence of the sub- 
leading terms neglected in the HTL-limit is considered. For the dispersion relations of the 
longitudinal and transverse gluons we find a decreasing gluon energy uj{p) with increasing 



g. This qualitatively resembles to the next-to-leading order behavior of the gluon mass |]15 
The deviations from the HTL result are larger than in the case of photons for the same 
extrapolated coupling strength (Fig. 8). 

IV. FERMION SELF ENERGY 

Here we consider the electron or quark self energy within the one-loop approximation. 
At zero temperature it is given by 

S(P) = zg'Cn J ^,DMYS{K)Y, (14) 

where the Casimir invariant of the fundamental representation is Cr = 4/3 for quarks 
(for electrons, g^C^ e^). We will first choose the Feynman gauge for the gauge boson 
propagator, 

D^^uiQ) = -9,.u-^ = -g^.u^B{Q) . (15) 

Combining (p!4D and ([T5|) , the finite temperature self energy is given by 

^ TY:I^A^{K)Ab{Q). (16) 

For evaluating this expression we utilize a general representation of the fermion self 
energy in the rest frame of the heat bath 

S(P) = -af - 670 (17) 

with the scalar functions 

a(P) = ^(tr[fS]-potr[7oS]) 

6(P) = ^ (p2tr[7oS] - potrirn) • (18) 



Evaluating the traces 

ri(p) = M7oS(P)], 

niP) = trirnp)] , (19) 

we are able to determine the functions a and b. Performing the sum over ko we get 

r,(P) = -89'C«/^ J^J, (20) 



where 



4b = I [1 + Mq) - nF(k)] ( ] + \ — 



1. 



+ [riBiq) + npik)] ( \ + \ ] \ . (21) 

^ ^ ^^\Po + k-q Po-k + qjj ^ ^ 

Again we neglect the vacuum contribution. After integrating analytically over the angle and 
numerically over k we end up with the final result for Ti. In Fig. 9 the real part of Ti is 
compared to the HTL result 

Tr^(P) = 2m2i In (22) 
^ P Po-P 



which contains the effective, thermal fermion mass 

mj^Cn^. (23) 



As for the boson self energies we observe the one-loop results approaching the HTL predic- 
tions for po,p <t^T and po ^ p (Fig. 9). 

Similar results can be found for T2, which follows from 

UP) = Sg'Cn / [P • ki^^ - Po F^%] , (24) 



where 



Akq [ \po + k + q PQ — k — q) 

+ [nB{q)+nF{k)]{ \ \ )l. (25) 



The comparison of the real part of T2 with the HTL result 

T^^^{P)=Am} (26) 

is shown in Fig. 10. 

In contrast to the HTL approximation, both one-loop fermion self energy functions de- 
velop an imaginary part. As shown in Fig. 11 and Fig. 12, these may be as large as the real 
parts. We note that the relevant combinations ±T2(P) — (±po — p)Ti{P) (see below) have 
non-positive imaginary parts. 

Next we consider the dispersion relations following from the resummed fermion propa- 
gator. It is convenient to decompose this propagator according to its helicity eigenstates 



with 



^ ~ 2D+(P) + 2D4P) ^^^^ 



D±{P) = -po ± p + 1 (±T2(P) - (±po - p)Ti(P)) . (28) 



The dispersion relations following from the roots of ( [28| ) are depicted in Fig. 13. The lower 
branch (plasmino) corresponding to D_(P) and displaying a minimum at finite momentum, 
which has a negative helicity to chirality ratio |jl3[, is absent in the vacuum. In Fig. 13 the 



complete one-loop dispersion relations (in Feynman gauge) are compared to the HTL ones. 
The differences increase with increasing coupling constant. However, even for g = "i the 
maximum deviation is of the order of 10%. 

Finally we investigate the gauge dependence of the one-loop fermion self energy. For this 
purpose we adopt the gauge boson propagator in a general covariant gauge 

D,,{K) = -g,,^ + ^K,K,j^ (29) 

with the gauge fixing parameter ^ (see e.g. [Q)- For treating the double pole in the second 
term we introduce a regularization parameter A ||16| and write 
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I _ d 1 



A=0 



d_ 



(30) 



A=0 



Here Kx indicates the four momentum (ko, kx) with \kx\ = kx = Vk"^ + \. Substituting the 



gauge dependent part of the boson propagator (^) into the self energy (jWj ) we obtain 



d^k 



3 Tj2mAB{Kx)Ap{Q) 

ko 



(31) 



A=0 



After some manipulations similar to the one for computing Ti and T2 we find the gauge 
dependent parts 6T1 and 6T2 of the complete one-loop expressions for Ti and T2 as a function 
of po and p (Fig. 14 and Fig. 15). The fact that by Ward-identities and power counting 
arguments they are of the same order as typical corrections of HTL by sub-leading 
terms is numerically found to hold also for hard momenta po,p ~ T. As expected, the gauge 
dependent contributions vanish for small po and p, where the one-loop self energy approaches 
the HTL result. 



V. CONCLUSIONS 

In this work we have calculated the complete one-loop self energies for photons, gluons, 
and fermions (electrons or quarks) and compared them to their HTL limits. In the case 
of the gauge independent one-loop photon self energy we found a good agreement between 
these two approaches. Especially for the dispersion relations both approximations agree to a 
surprisingly high accuracy. As far as the real parts of the gluon and fermion self energies are 
concerned, the difference between the both approaches is more distinct. For momenta near 
the light cone, however, both approaches coincide even for hard momenta violating the HTL 
assumption pQ,p <^T. This means that for the HTL approximation being valid we need not 
to require that both pq and p are soft, but assuming a soft P"^ = Pq — p^ is sufficient. This 
also justifies the use of the HTL resummed quark propagator for soft four- momenta in the 
calculation of the production rate of hard photons from the QGP ||T^. For the imaginary 
parts there are somewhat larger differences between the complete one-loop and the HTL 
calculation. In particular the one-loop fermion self energy shows a significant imaginary 
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part above the light cone which may be of the same order as the real part. (The same is 
expected for the (complete) gluon self energy beyond leading HTL order.) Finally we have 
studied the gauge dependence of the one-loop fermion self energy in the class of covariant 
gauges. 

The results found above indicate that the HTL approximation, based on the weak cou- 
pling limit assumption (7 ^ 1 or equivalently on the high temperature assumption T ^ Po,P, 
is qualitatively correct for a large momentum regime even for large values of (? > 1 or tem- 
peratures of the order p. This gives us some confidence in the usefulness of the effective 
perturbation theory (Braaten-Pisarski method), which is based on the HTL resummation, 
for realistic values of the coupling constant g = 1.5 - 2.5 at temperatures T = 200 - 300 
MeV that are accessible in relativistic heavy ion collisions. 

Furthermore, the use of HTL effective masses for describing the equation of state of a 
gluon gas m or strange quark matter p[ appears to be justified from the point of view 
that the deviation of the HTL masses from the complete one-loop effective masses can be 
neglected in those models. 
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FIG. 1. Real part of the longitudinal photon self energy within the complete one-loop approx- 
imation (solid line) and the HTL limit (dashed line) for momenta p/T = 0.5, 1, 1.5. 
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FIG. 3. Dispersion relation of longitudinal and transverse photons (lower and upper branches, 
respectively) within the complete one-loop approximation for e = 3 (solid line) and the HTL limit 
(dashed line). 
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PolT 

FIG. 5. Real part of the transverse gluon self energy within the complete one-loop approxima- 
tion in Feynman gauge (solid line) and the HTL limit (dashed hne) for p/T = 0.5, 1, 1.5. 
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FIG. 6. Imaginary part of the longitudinal gluon self energy within the complete one-loop 
approximation in Feynman gauge (sohd hne) and the HTL limit (dashed line) for p/T = 0.5, 1, 1.5. 
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FIG. 7. Imaginary part of the transverse gluon self energy within the complete one-loop ap- 
proximation in Feynman gauge (solid line) and the HTL limit (dashed line) forp/r = 0.5,l,1.5. 
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FIG. 8. Dispersion relation of the longitudinal and transverse gluon modes (lower and upper 
branches, respectively) in one-loop approximation in Feynman gauge for 5 = 3 (solid line) and 
HTL hmit (dashed). 
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FIG. 10. Real part of T2 within the complete one-loop approximation in Feynman gauge (solid 
line) and the HTL limit (dashed line) for p/T = 0.5, 1, 1.5. 
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FIG. 11. Imaginary part of T\ within the complete one-loop approximation in Feynman gauge 
(solid line) compared to HTL result (dashed line) for -pjT = 0.5, 1, 1.5. 
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FIG. 12. Imaginary part of T2 within the complete one-loop approximation in Feynman gauge 
(solid line) for p/T = 0.5, 1, 1.5. Tj^"^^ has no imaginary part. 
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FIG. 13. Quark dispersion relations within the complete one-loop approximation in Feynman 
gauge for 5 = 3 (solid line) and the HTL limit (dashed line). The upper branch is the particle 
exitation, the lower one is the plasmino. 
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FIG. 15. Gauge dependent contribution of the the real part of T2 {p/T = 0.5, 1, 1.5). 
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